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Abstract. We discuss the details of calculating hadron properties from the OPE for correlators 
of quark currents in QCD, which constitutes the basis of the method of QCD sum rules. 
The main emphasis is laid on gaining control over the systematic uncertainties of the hadron 
parameters obtained within this method. We start with examples from quantum mechanics, 
where bound-state properties may be calculated independently in two ways: exactly, by solving 
the Schrodinger equation, and approximately, by the method of sum rules. Knowing the exact 
solution allows us to control each step of the sum-rule extraction procedure. On the basis of 
this analysis, we formulate several improvements of the method of sum rules. We then apply 
these modifications to the analysis of the decay constants of heavy mesons. 



1. Introduction 

The method of dispersive sum rules yij is one of the widely used methods for obtaining properties 
of ground-state hadrons in QCD. The method involves two steps: (i) one calculates the relevant 
correlator in QCD at relatively small values of the Eucledian time; (ii) one applies numerical 
procedures suggested by quark-hadron duality in order to isolate the ground-state contribution 
from this correlator. These numerical procedures cannot determine a single value of the ground- 
state parameter but should provide the band of values containing the true hadron parameter. 
This band is a systematic, or intrinsic, uncertainty of the method of sum rules. 

An unbiased judgement of the reliability of the extraction procedures adopted in the method 
of sum rules may be acquired by applying these procedures to problems where the ground-state 
parameters may be found independently and exactly as soon as the parameters of theory are 
fixed. Presently, only quantum-mechanical potential models provide such a possibility: (i) the 
bound-state parameters (masses, wave functions, form factors) are known precisely from the 
Schrodinger equation; (ii) direct analogues of the QCD correlators may be calculated exactly. 

Making use of these models, we studied the extraction of ground-state parameters from 
different types of correlators: namely, the ground-state decay constant from two-point vacuum- 
to- vacuum correlator [2] , the form factor from three-point vacuum-to- vacuum correlator pj , and 
the form factor from vacuum-to-hadron correlator [3] . We have demonstrated that the standard 
procedures adopted in the method of sum rules not always work properly: the true value of the 
bound-state parameter was shown to lie outside the band obtained according to the standard 
criteria. These results gave us a solid ground to claim that also in QCD the actual accuracy 
of the method may be worse than the accuracy expected on the basis of applying the standard 
criteria. 



We realized that the main origin of these problems of the method originate from an over- 
simplified model for hadron continuum which is described as a perturbative contribution 
above a constant Borel-parameter independent effective continuum threshold. We introduced 
the notion of the exact effective continuum threshold, which corresponds to the true bound- 
state parameters: in potential models the true hadron parameters — decay constant and form 
factor — are known and the exact effective continuum thresholds for different correlators may 
be calculated. We have demonstrated that the exact effective continuum threshold (i) is not a 
universal quantity and depends on the correlator considered (i.e., it is in general different for 
two-point and three-point vacuum-to- vacuum correlators), (ii) depends on the Borel parameter 
and, for the form-factor case, also on the momentum transfer [31lll[5]. 

In recent publications [6j we proposed a new algorithm for extracting the parameters of the 
ground state. The idea formulated in these papers is to relax the standard assumption of a 
Borel-parameter independent effective continuum threshold and to allow for a Borel-parameter 
dependent quantity. This lecture explains the details of this procedure and its application to 
decay constants of heavy mesons. 

2. OPE and sum rule in quantum- mechanical potential model 

Let us start with quantum-mechanical potential model described by Hamiltonian 

A;2 a 
H = Ho + V{r); Ho = —, V{r) = V,ondr) - (1) 
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Polarization operator Il{E) and its Borel transform n(T) [E — )• T, 1/{H — E) ^ exp{—HT), T 
the Borel parameter] are defined via the full Green function G{E) = 1/{H — E) [7j: 

n(^) = {rj = ^\jPzr^\n = 0), n(r) = (r) = o| eM-HT)\n = o). (2) 

The expansion of G{E) and !!(£') in powers of the interaction is obtained with the help of the 
Lippmann-Schwinger equation 

G{E) = Gq{E) - Go{E)VGo{E) + Go{E)VGo{E)VGo{E) + ..., (3) 

where Go{E) = 1/{Hq — E). Since the interaction contains now two parts, V"conf('') and ^, the 
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expansion ([3]) is a double expansion in powers of V^onf and a. E.g., for the case Konf = 
one easily obtains the corresponding double expansion in powers of a and ojT: 



noPE(r) = npert(T)+npower(T), ^P^rt{T)= i^j 
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(4) 



One can see here a "perturbative contribution" (i.e. not containing the confining potential), 
and "power corrections" given in terms of the confining potential (including also mixed 
terms containing contributions from both Coulomb and confining potentials). A perturbative 
contribution may be written in the form of spectral representation [8j yielding 



dze~^'^ Ppert{z) + npower(7'), Ppert = (^) ^Vvt"^ V^^Jima + 



(5) 



The "physical" representation for n(T) — in the basis of hadron eigenstates — reads: 

oo 

nphy.(T) = {ff = 0| exp(-FT)|f, = 0) = ^ RneM-EnT), i?„ = |^'„(f = 0)\\ (6) 

n=0 

Sum rule is the expression of the fact that the correlator may be calculated in two ways — using 
the basis of quark states (OPE) or confined bound states — leading to the same result: 

noPE(r) = nphys(r). (7) 

In order to isolate the ground-state contribution one needs the information about the excited 
states. A standard Ansatz for the hadron spectral density has the form [Ij 

Pphys{z) = Rgd{z - Eg) + e{z - Zc{i)ppcTtiz). (8) 

It assumes that the contribution of the excited states may be described by contributions of 
diagrams of perturbation theory above some effective continuum threshold Zes- This effective 
continuum threhsold (different from the physical continuum threhsold which is determined by 
hadron masses) is an additional parameter of the method of sum rules. Using Eq. ([5]) yields 

^cff 

Rge-^<'^ = J dze-'^ppertiz) + npower(T) = Ud^^T, ^eff)- (9) 


As soon as one knows Zcs, one immediately obtains estimates for Rg and Eg, R^naiiT, Zef[) and 
EduaiiT, Zes) = — (ir log n(T, Zeff). These however depend on unphysical parameters T and ZeS- 

Eq. dS]) is motivated by quark-hadron duality which claims that far above the threshold the 
hadron spectral density is well described by diagrams of perturbation theory. However, near 
the physical threshold — and this very region turns out to be essentail for the calculation of 
ground-state properties — the duality relation is violated. 

The advantage of quantum mechanics is that the exact Eg and Rg may be obtained by solving 
Schrodinger equation [S]. We can then calculate z^g from 

Rge-^^^ = y dze-^^ppert(2)+npowcr(T). (10) 



The obtained "exact threshold" 2;efj(r) is a slightly rising function of T [2]- 



3. OPE and sum rule in QCD 

In QCD, the sum rule for the decay constants of heavy mesons has a similar form jlUj 

/q^QE"^^^ = J e"''^Ppert(s, a, 'mQ,p) ds + npower(-r, mQ,p) = ndual(T", •Seflf) (11) 
(mQ+r?iu)2 

In order to extract the decay constant one should fix the effective continuum threshold Seff which 
should be a function of r; otherwise the r-dependences of the l.h.s. and the r.h.s. of ()lip do 
not match each other. The exact Seg corresponding to the exact hadron decay constant and 
mass on the l.h.s. is of course not known. The extraction of hadron parameters from the sum 
rule consists therefore in attempting (i) to find a good approximation to the exact threshold 
and (ii) to control the accuracy of this approximation. For further use, we define the dual decay 
constant /dual ^ind the dual invariant mass Mjuai by relations 

/d'ual(T) = MQ^e^'Q^nduai(r,Sefr(T)), Mi,i(r) = - A log nduai(r, Sefr(T)). (12) 



4. Extraction of the decay constant: QCD vs. potential model jllj 

If the mass of the ground state is known, the deviation of the dual mass from the actual mass of 
the ground state gives an indication of the contamination of the dual correlator by excited states. 
Assuming a specific functional form of the effective threshold and requiring the least deviation 
of the dual mass ()12p from the known ground-state mass in the r-window leads to a variational 
solution for the effective threshold. As soon as the latter has been fixed, one calculates the decay 
constant from (|12|) . Our algorithm for extracting ground-state parameters reads: 



(i) Consider a set of Polynomial r-dependent Ansaetze for Ses'- s^^{t) = Sq'' +Si''t+S2'''t''^+- ■ ■ ■ 

(ii) Calculate the dual mass for these r-dependent thresholds and minimize the squared difference 
between the "dual" mass M^^^^ and the known value in the r-window. This gives the 

(n) 

parameters of the effective continuum thresholds s] . 

(iii) Making use of the obtained thresholds, calculate the decay constant. 

(iv) Take the band of values provided by the results corresponding to linear, quadratic, and 
cubic effective thresholds as the characteristic of the intrinsic uncertainty of the extraction 
procedure. 

Figure [1] shows the application of our algorithm in quantum mechanics and in QCD. In both 
cases a very similar pattern emerges; the reproduction of the dual mass considerably improves 
for the r-dependent quantities indicating that the dual correlator with a r-dependent threshold 
isolates the contribution of the ground state much better than the dual correlator with a standard 
r-independent threshold. As a consequence, the accuracy of the extracted hadron observable 
improves considerably. 
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Figure 1. The outcome of our algorithm in potential model (left) and in QCD (right). 



5. Decay constants of D and Ds mesons [12j 

The application of our extraction procedures leads to the following results: 

fD = 206.2 ± 7.3(oPE) ± 5.1(3y,t)MeV, fo, = 245.3 ± 15.7(ope) ± 4.5(,y,t) MeV. (13) 

The r-dependent threshold is a crucial ingredient for a successful extraction of the decay constant 
from the sum rule (Fig. [2]). Obviously, the standard r-independent approximation leads to much 
lower value for /£> which lies rather far from the data and the lattice results. 
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Figure 2. Our results for and f^s [12] , 

6. Conclusions 

The effective continuum tiireshold s^s is an important ingredient of the method of dispersive sum 
rules which determines to a large extent the numerical values of the extracted hadron parameter. 
Finding a criterion for fixing Sgfr poses a problem in the method of sum rules. 

• Seff depends on the external kinematical variables (e.g., momentum transfer in sum rules for 3- 
point correlators and light-cone sum rules) and ' 'unphysical" parameters (renormalization scale 
fj,, Borel parameter r). Borel-parameter r-dependence of Ses emerges naturally when trying to 
make quark-hadron duality more accurate. 

• We proposed a new algorithm for fixing r-dependent Ses, for those problems where the bound- 
state mass is known. We have tested that our algorithm leads to more accurate values of 
ground-state parameters than the "standard" algorithms used in the context of dispersive sum 
rules before. Moreover, our algorithm allows one to probe "intrinsic" uncertainties related to 
the limited accuracy of the extraction procedure in the method of QCD sum rules. 
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